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The transport properties of Nb/CuMn multilayers with 
a regular array of electron beam lithography obtained an- 
tidots have been measured at different temperatures in the 
presence of external perpendicular magnetic fields. Hystereti- 
cal I ~ V characteristics have been observed which disap- 
pear when approaching the upper critical magnetic field curve 
Hc2{T). Comparing these data with other results (Arrhenius 
plots of resistive transition curves, logV-logI characteristics) 
we have been able to relate the onset of the hysteresis to the 
presence of an irreversibility line. We discuss several possible 
mechanisms to clarify the nature of this line. Among them the 
most plausible seems to be the vortex melting mainly induced 
by quantum fluctuations. 

PACS: 74.60.Ge, 74.60.Jg. 



I. INTRODUCTION 

The task of increasing the critical current density Jc 
in superconducting materials has always been a widely 
studied subject gaining even more interest after the 
discovery of high temperature superconductors (HTS) 
These studies are strictly related to the knowledge 
of the flux line pinning mechanism and to the reduction 
of the vortex mobility. Introducing artificial defects in 
superconducting materials as, for example, non super- 
conducting distributed phases ||,^, columnar tracks of 
amorphous material obtained by high energy ions H or 
geometrical constrictions such as channels or dots 
is a very useful tool for a better understanding of the 
vortex dynamics and for obtaining higher Jc values. The 
recent developments of the submicrometer electron beam 
lithographic techniques have made possible to reduce the 
typical size of these geometrical constrictions to values 
much smaller than the period of the vortex lattice, and 
comparable to those of typical superconducting coher- 
ence lengths P,pHlO[. Many experiments, performed on 
systems with such a regular array of defects and 
several numerical simulations , have been focused 

on the study of the vortex properties at low magnetic 
fields close to the matching fields if„ = rip^o (here Up is 
the pins concentration and $o is the flux quantum). 



Other studies, performed at higher magnetic fields, 
have been related to the analysis of the vortex lattice 
shear stress in superconducting layered systems with the 
presence of artificially obtained weak pinning channels 
embedded in a strong pinning environment |^,|^,|l^,|l6| . 
Due to the possibility of varying in a controlled way many 
different parameters such as the density, the dimension- 
ality and the nature of the pinning centers , the 
study of the vortex dynamics in artificially layered con- 
ventional superconducors with the simultaneous presence 
of a regular array of pinning centers perpendicular to the 
layers is of great interest. Moreover, the natural layered 
structure of HTS allows to use these artificial systems as 
a model to help in discriminating among intrinsic and 
dimensional effects in the transport properties of HTS 
compounds Q. 

In this paper we report on current- voltage {I—V) char- 
acteristic and resistance versus temperature, R(T,H), 
measurements, in perpendicular external magnetic fields 
H , performed on two different series of Nb/CuMn mul- 
tilayers with a square array of antidots. The choice of a 
superconducting (Nb)/spin glass (CuMn) layered system 
could be interesting particularly in view of the use as a 
model of HTS compounds. The two series have antidots 
with the same diameter, I? « 1 ^m, and different lat- 
tice distances between the antidots: d « 2 /im in one 
series and d « 1.6 /im in the other. The experiments 
have been performed on different samples having different 
anisotropics for each series. The I—V curves measured at 
different temperatures and at different values of H have 
shown, in regions of the H — T phase diagram depending 
on the anisotropy of the system, a hysteretic behavior 
with sudden voltage jumps, which disappers approach- 
ing the critical field Hc2{T) curve. From the analysis 
of the curvature of the logarithmic I—V characteristics 
and from the study of the shape of the Arrhennius plots 
of the resistive transition curves, we have been able to 
relate the disappearing of such a hysteretical behavior 
to the presence of an irreversibility line (IL). We have 
discussed different possible mechanisms responsible for 
this IL. Among them the most plausible for our samples 
seems to be the vortex melting mainly induced by quan- 
tum fluctuations ||l9| . 
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II. EXPERIMENTAL 

Regular arrays of antidots have been fabricated using 
Electron Beam Lithography to pattern the resist on a 
2 inch Si (100) wafer. The Nb and CuMn layers were 
deposited by a dual-source magnetically enhanced dc tri- 
ode sputtering system with a movable substrate holder 
onto the patterned resist and the final structures were 
obtained by the lift-off technique. Single antidots have 
a circular geometry and the antidot array is arranged in 
a square lattice configuration, see figure |l|. The total 
area covered by the array is a square of 200 x 200 //to^ 
with four separate contact pads connected to the vertices. 
Eight replicas of this structure are present on the same Si 
wafer to allow the fabrication of a series of multilayered 
samples with the same Nb thickness and variable CuMn 
thicknesses in only one deposition run [Q. The resist 
used is UV III from Shipley; UV III is a chemically am- 
plified photoresist for the deep-UV range, but is widely 
used as an e-beam resist because it provides a good trade- 
off between reasonable sensitivity and high resolution. In 
our case it has been used as a positive resist, that is, the 
resist is retained where unexposed. A resist film thick- 
ness of 5800 A has been obtained spinning the wafer at 
1800 rpm. The resist has been exposed using a Leica 
Cambridge EBMF 10 system operated at 50 kV. Several 
tests have been carried out in order to achieve structure 
profiles most suitable for lift-off. In particular the de- 
sired profile is that showing a moderate undercut. Also, 
developing times and post-exposure treatment have been 
optimized for profile improvement. After developing and 
post-baking at 130 °C the samples underwent RIE in O2 
for 30 s at 25 W rf power and 14 Pa oxygen pressure, 
in order to completely remove residual resist in the ex- 
posed areas. This treatment lowered the resist thickness 
to about 5000 A. Fig.l shows a Scanning Electron Micro- 
scope (SEM) image of the typical result of the fabrication 
process: in this case the antidots nominal diameter is 1 
jjLm and the nominal period of the structure is 2 jjun. 
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Table I. Some of the relevant sample features, pat is 
the resistivity at T=10 K. See the text for the meaning 
of the other quantities. 



The number of bilayers of Nb and CuMn is always 
equal to six. The first layer is CuMn, the last one is Nb. 
The Nb nominal thickness, dNb, is 250 A for all the sam- 
ples in both the series. The CuMn thickness, dcuMn, has 
been varied from 7 A to 25 A in one series (d « 2/Ltm) and 
from 4 A to 20 A in the other {d « 1.6^m). The Mn per- 
centage is always equal to 2.7. For reference a Nb/CuMn 
multilayer with dc'uMn=13 A (Mn %=2.7) without an- 
tidots, but with the same configuration (200x200 fim^ 
square and four pads connected to the vertices), has 
been fabricated. The samples parameters are summa- 
rized in Table 1. I — V characteristics have been registred 
at T < 4.2 K using a dc pulse technique. The tempera- 
ture stabilization, during the acquisition of curves in the 
helium bath, was better than 10^^ K. The magnetic field 
was obtained by a superconducting Nb-Ti solenoid. From 
the measured temperature dependencies of the perpen- 
dicular and parallel upper critical field Hc2 (obtained at 
half of the resistive transitions R{T, H)) we deduced the 
values for the parallel and the perpendicular coherence 
length at zero temperature, C||(0) and ^±(0) respectively, 
and then the value of the anisotropic Ginzburg-Landau 
mass ratio 70 = C_l(0)/^|| (0) ||2l|. The different values 
of the critical temperatures for the samples of the two 
series having similar values of dNb and dcuMn are prob- 
ably related to the different Nb quality obtained in the 
two deposition runs. 

Figure |^ shows I — V curves for the sample NCMF 
at r=2.60 K for different applied magnetic fields in the 
range 0.03< H <0.70 Tesla. At low magnetic fields 
the curves present hysteresis, not due to thermal effects, 
when registered both in forward and backward directions, 
i.e. increasing and decreasing the current ||l^. Such hys- 
teresis becomes smaller when H is increased disappear- 
ing completely, in the limit of our experimental accuracy, 
when approaching Hc2- These features have been repeat- 
edly obtained for the same sample in different measure- 
ments and are typical for all the samples of the two series 
with antidots. In figure ^3 the I — V characteristics for 
the sample NCMD are plotted in the temperature range 
2.3 K< T <2.91 K for ^ioH^ 0.03 T. In this case the hys- 
teresis disappears approaching T^. The I — V curves for 
the sample NCMDA without antidots are always smooth 
and parabolic-like, typical of a type II superconductor. 

In figure |^, we plot the I — V curves in double loga- 
rithmic scale for the sample RNCMC at r=2.96 K for 
different values of the magnetic field. A change in the 
curvature of the logV-logI curves clearly occurs at low 
voltages; this is usually related to the presence of an IL 
in the H — T phase diagram [ p2[ , given by the points 
at which the logV-logI curve is linear. In all the sam- 
ples measured, these points are always very close to the 
points at which the hysteresis in the I — V curves disap- 
pears. As an example, in the figures ^(a,b,c) are shown 
the H-T phase diagrams for the samples NCMB, NCMC 
and NCMH, respectively. Circles distinguish the two dif- 
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FIG. 1. Scanning Electron Microscope picture of a 
Nb/CuMn sample with a square lattice of antidots {D ~ 1 
fim and d ~ 2 ^m). 
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FIG. 3. Logarithmic Current- Voltage characteristics for 
the sample RNCMC at T = 2.96 K for magnetic fields, in- 
creasing counterclockwise, of 0.20; 0.40; 0.50; 0.80; 0.90; 1.10; 
1.30; 1.50; 1.60; 1.70; 1.75; 1.80; 1.85; 1.95; 2.10 T. The num- 
bers indicate some values of the magnetic field. The melting 
field is estimated to be equal to 1.75 T. 
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FIG. 2. (a) Current- Voltage characteristics for the multi- 
layer NCMF at T = 2.60 K for different magnetic field values. 
The increasing values are: 0.03; 0.06; 0.15; 0.55; 0.70 T. (b) 
Current- Voltage characteristics for the multilayer NCMD at 
H = 0.03 T for different temperatures. The increasing values 
are: 2.30; 2.34; 2.60; 2.80; 2.91 K. 



ferent regimes of the vortex dynamics as determined from 
the disappearance of the hysteresis in the 1~Y curves; up 
triangles indicate points where change in the curvature 
of the logV-logI occurs; squares correspond to the values 
of the perpendicular critical magnetic field. It is then 
evident that the change in the hysteretical behavior of the 
1—V characteristics, is related to the IL as defined by the 
change in the logV-logI curvature. It is also interesting 
to note that the position of the IL in the H — T plane 
moves a'way from the 77^2 line as the anisotropy of the 
samples increases. 

Another "way to determine the presence of an IL in the 
H — T phase diagram of a superconductor is the study of 
the Arrhenius plot of the resistance versus the temper- 
ature curves [|l8| . In figure || the Arrhenius plot of the 
transiton curves, recorded using a bias current of 2 niA, of 
the sample RNCMA is sho-wn at different perpendicular 
magnetic fields. A 'well defined field dependent temper- 
ature T* separates t'wo zones 'with very different activa- 
tion energies. In particular, at T < T* a sudden increase 
in the Arrhnius slope signals a transition in the trans- 
port properties of the sample ■which can be related to the 
presence of the IL. In figure ^ it is sho'wn the measured 
H — T phase diagram for the sample RNCMC. The solid 
squares correspond to the perpendicular magnetic fields; 
the open squares correspond to the points in the H — T 
plane at -which the onset of hysteresis in the I —V curves 
takes place, the open diamonds are defined taking into 
account the change of curvature of the logl-logV curves, 
and the open circles and up triangles are the points at 
-which the slope in the Arrhenius plot of the transition 
curves, taken respectively at a bias currents of 2 mA and 
6 mA, change as vie'wed in figure ^. It is evident that 
also in this case all the points (open symbols) fall again 
on the same curve, 'which can be identified as an IL. 
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III. DISCUSSION 

In all the measurements performed the bias current 
was applied in the plane of the film and perpendicular 
to the direction of the external magnetic field, see fig- 
ure 1^. In this configuration, the Lorentz force acting on 
the vortices, tends to move them along the channels in 
between the antidots. On the other hand, in the nar- 
rower zones between adjacent antidots, the current den- 
sity is much higher than in the channels, so that we lo- 
cally have weaker pinning centers in these parts of the 
sample. Therefore, we can look at our superconducting 
layered system as made of alternating zones of strong 
pinning (the channels along the direction of action of the 
Lorentz force in figure and weak pinning (the narrower 
zones between adjacent antidots], similarly to other cases 
reported in the literature [|6P Jl5| , |l6|] . The value of the 
matching field in both the series is very small, w 5 
Oe and also d ^ £,{T), where £,{T) is the temperature 
dependent coherence length which in our multilayers has 
typical values of about 100 A. Therefore, in the supercon- 
ducting state it is possible to have a vortex lattice inside 
the channels between the antidots. The pinning of these 
interstitial vortices is determined mostly by the intrinsic 
properties of the superconducting materials ■ 

To start the analysis of our data we have first to de- 
termine the dimensionality of the vortex lattice in our 
samples. At magnetic fields lower than a critical value 
Her ~ 4$o/7o'S^, where s is the interlayer distance be- 
tween superconducting layers, the vortices in adjacent 
layers are strongly coupled and vortex lines behave as 
three domensional (3D). For our samples Her ~ 10^ T, 
well above any applied fields, and we can exclude the ap- 
pearance of decoupling of the vortex lines in our layered 
systems Q. 

On the other hand, when the shear modulus Cge of a 
vortex lattice is much lower than the tilt modulus C44, the 
system can behave bidimensionally. In principle, thermal 
fiuctuations could cause tilt deformations in the vortex 
line. The critical thickness value dcr of a film beyond 
which thermal fluctuation induced tilt deformations be- 
come relevant is given by |^ 



4.4^11 



(1) 



where h = H/Hc2- In multilayers, however, the dcr 
value is reduced by a factor 7q which for our samples 
gives, in the range of measured temperatures and mag- 
netic fields, dj."^""* « 200 300 A. Typical thickness for 
our samples are in the range 1500 1700 A. This means 
that the vortex lattice in our multilayers is in a strong 
3D regime. 
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FIG. 4. H-T phase diagram for the sample (a) NCMB, (b) 
NCMC and (c) NCMH. Squares correspond to Hc2±, circles 
correspond to the disappearence of hysteresis in the I — V 
curves, up triangles correspond to the change of the curva- 
ture in the logV-logI curves. The solid lines are calculated 
according to the equation 3. 
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FIG. 5. Arrhenius plot of the transition curves for the sam- 
ple RNCMA in perpendicular magnetic fields. The increasing 
magnetic field velues are: 1.2; 1.4; 1.6; 1.8; 2.0; 2.2; 2.4; 2.6; 
2.8 T. The numbers indicate some values of the magnetic field. 
The curves have been recorded using a bias current Ib= 2 mA. 
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FIG. 6. H-T phase diagram for the sample RNCMC. Solid 
squares correspond to Hc2±, open squares correspond to the 
disappearence of hysteresis in the I ~ V curves, diamonds 
correspond to the change of the curvature in the logV-logI 
curves, circles (up triangles) correspond to the T*{H) values 
extracted from the Arrhenius plot using a bias current of 2 
mA (6 mA). The solid line has been calculated according to 
the equation 3. 

One of the proposed interpretation of the nature of 
the IL relies upon a depinning mechanism in which a 
crossover from flux creep to flux flow occurs ||26|] . In our 
samples we never obtain linearity of the I —V curves nei- 
ther at high voltage where a uniform flux flow should be 
present, nor at small currents, where thermally assisted 
flux flow (TAFF) d?) should take place. If we flt our data 
with the relation V ^ I" the fltting exponent is always 
very high (a > 10). Therefore we exclude the possibil- 
ity of flux creep-flux flow crossover in the presence of a 
pinning strength distribution which could also be respon- 
sible for the curvature change of the I — V curves when 



a and 0b are very 



plotted in double logarithmic scale pg|. 

The I — V curves shown in figure 
similar to those obtained as a result of numerical simu- 
lations for a superconductor with periodic pinning close 
to the matching field |13 1^. Considering that we are 
very far from Hn we cannot apply the results of Ref. |l^] 
to explain our I — V curves. Nevertheless, the main ob- 
served features indicate the presence of a region in the 
H — T plane where plastic vortex motion probably takes 
place. Below some crossover value Hpi, vortices experi- 
ence plastic motion which usually reveals in hysteretical 
curves . Above Hpi , I — V curves are smooth, hys- 
teresis vanish and the vortex motion start to be a fiow 
of vortex liquid |^9|] . The vortex melting scenario along 
with the simultaneous presence of weak and strong pin- 
ning channels, can also explain the observed Arrhenius 
plots. In fact, immediately below J'°"«'=* (defined as the 
temperature where the electrical resistance R is 0.9 of its 
value in the normal state), the solidification of the vor- 
tices in the strong pinning channels determines the high 
values of the initial slope in the plots. At lower tempera- 
tures, the dissipation in the system is mainly due to the 
vortices in the weak pinning regions, and this results in 
a lower value of the activation energies. When also these 
vortices experience the transition from liquid to solid at 
T = r* , the slope in the plots increases again ||l6| . 

Melting in the vortex lattice can be induced by thermal 
fiuctuations ||3^. The melting temperature at which the 
vortex lattice goes from a solid phase to a liquid one, can 
be obtained by using the 3D thermal melting criterion 
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where cl is the Lindemann number, Tc is the supercon- 
ducting transition temperature (defined in our case at 
the point where the electrical resistance R of the sam- 
ple becomes less than 10"'' ft), t„i — Tm/Tc is the 
reduced melting temperature, and Gi is the Ginzburg 
number which determines the contribution of the ther- 
mal fluctuations to the vortex melting given by Gi — 
(l/2)(27r^oA:BT,A||(0)7o/$^^||(0)), where \\\ is the in- 
plane penetration depth which, for all our samples, has 
been assumed to be equal to 1500 A ||3l[ . Melting usu- 
ally occurs when cl ~ 0.1 -^0.3 ]19(|. When we try 
to flt the IL observed in our samples with Eq.2 we get 
Cl ^ 10^^. This extremely low value for the Lindemann 
number makes unreasonable to consider the IL as due to 
3D thermal melting. 

However, as first pointed out by Blatter and Ivlev |l^] 
in moderately anisotropic superconductors at low tem- 
peratures one could not exclude the contribution of quan- 
tum fluctuations to the melting. In this case the total 
fluctuation displacement of vortex line is < u >^=< 
.2 n, \2 |--|^g average dis- 



>t/i + < " >g; where a/< u > 
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placement due to thermal fluctuations and y < u >^ is 

the average displacement due to quantum fluctuations. 
< u diminishes with the temperature, while < u >^ 
does not depend on the temperature and at low values 



ad H,2{T) 
k2 r(1.25 - 0.25i) 



(1 - 0.58/i„ - 0.29/i2j(l - h„,f = 1 

(4) 



of T one could expect < u >t^<< u The ampli- we do not obtain any agreement with the data. Here 

A$o(1.07)2/327r^ofcs, « = A||/^|| and d is the thick- 



tude of < M >g depends on the ratio Q* j^JTri where 
Q* — e^pisi/fLS, with U the Planck constant and e the 
elementary charge. If Q* l^fGl ^ 1, the contribution 
of quantum fluctuations is crucial. For the samples dis- 
cussed here we always get Q*/\/Gi > 30 which justifies 
the possibility of an important contribution coming from 
quantum fluctuations In this case the melting line 
is given by [pi 
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1 + (i + 4QeJ- 



(3) 



where hm = Hm/Hc2, t = T/Tc is the reduced temper- 
ature, e = 7rc|(t-i - l)/VU'i, Q = Q*nT/'K^/Gl with 

a cut-off frequency generally of the order of the De- 
bye frequency and r an effective electronic relaxation 
time [fi/kBTc ~ t) [Q. As we already showed be- 
fore [Q, the values of VL and r in Nb/CuMn are in the 
range (2^3) x lO^^s"^ and (1-^5) x lO^^^s respectively. 
Therefore, for all the samples studied we have assumed 
n = 3 X lO^^s"^ and t = 5 x 10~^^ s. In this way we 
reduce the number of free fit parameters in Eq.3 to only 
one, namely the Lindemann number cl- 

The solid line in Fig. 6 has been calculated according to 
Eq.3 using for the Lindemann number a value of cl=0.23. 
Also the solid lines in the Figures ^b and ^ have been 
calculated according to Eq.3. For the sample NCMC, 
Fig. 4b, we obtain good agreement with the experimental 
data for cl = 0.19, while for the sample NCMH, Fig. 4c, 
the solid fine is obtained taking cl = 0.09. The agree- 
ment between experimental data and theoretical curves 
is very good for all the samples studied. The Cl val- 
ues, as shown in Table I, become smaller with increasing 
anisotropy reaching in the case of sample NCMH a value 
slightly below 0.1 When increasing the anisotropy 

of the system, the coupling between adjacent supercon- 
ducting layers reduces and vortex lines become more soft. 
The influence of the thermal fluctuations on the vortex 
dynamics is strongly dependent on this coupling. Qual- 
itatively, therefore, the reduction of the cl values with 
increasing anisotropy could be related to the change in 
the topology of the vortex system. 

On the other hand, the softening of the vortex system 
could also determine a situation in which thermal fluctu- 
ations are able to cause tilt deformations. Anyway, if we 
try to fit the experimental points in figures ^ and ^ using 
the 2D pure thermal melting curve |3q| 



ness of the sample. A is a renormalization factor of the 
shear modulus Cge due to non linear lattice vibrations 
and vortex lattice defects and is ^ « 0.64 |36|. 
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FIG. 7. Values of the activation energy in the zone 2 for 
a not-perforated sample (solid squares) and values of the ac- 
tivation energy in the zone 1 for a perforated sample (open 
squares). Insert: schematic of a typical Arrhenius plot ob- 
served in a perforated and a not-perforated Nb/CuMn multi- 
layer together with the identification of the zones 1 and 2. 

We want to point out that the quantum melting theory 
has been succesfuUy applied to describe the vortex be- 
havior also in not-perforated Nb/CuMn multilayers [^ . 
In that case the melting line was determined by analyz- 
ing in Arrhenius fashion the measured R{T) curves in 
perpendicular magnetic fields. The shapes of the Ar- 
rhenius plots, see for example figure |l| in ref. p^ , were 
very similar to those observed in the case of perforated 
samples, suggesting the presence of two types of pinning 
centers also in the case of non perforated samples. In not 
perforated samples, edge pinning could be relevant and 
obviuosly stronger than intrinsic pinning j3^] . Therefore 
one could interpret the shape of the Arrhenius plots in 
not perforated Nb/CuMn multilayers as due to vortex 
transition from liquid to solid first of the vortices at the 
edges and then, at lower temperatures, of the vortices 
intrinsically pinned in the inner part of the samples. If 
this interpretation is correct, the slopes of the Arrhenius 
plots measured in not perforated samples in the zone 2 
(see insert in figure |^) should be very close to the slope 
measured in the zone 1 in perforated samples. In fact, 
in both cases, this slope should be related to the acti- 
vation energy of vortices intrinsically pinned inside the 
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system. In figure |7| the solid points refer to the values of 
the activation energy in a typical not-perforated sample 
measured in zone 2, while the open symbols refer to the 
values of the activation energy in a perforated sample 
(RNCMA) measured in zone 1. The two samples have 
been chosen to have similar critical temperatures. Also 
the R(T, H) curves have been takin using a similar value 
for the bias current density Jf,. The quite nice agreement 
between the two sets of data supports our idea that, in 
both cases, they are a measure of the intrinsic pinning in 
the material. 

The presence of the antidot array in the multilayers 
makes more easily measurable the vortex melting. In 
fact, in the case of perforated samples one is able to de- 
tect the change in the slope of the Arrhenius plots using 
low values of the bias current (~ 100/iA) while for not 
perforated samples bias currents of ~ 1mA are needed 
to observe the same effect (the IL in a superconductor 
does not depend on the value of the bias current). This 
is consistent with the idea that in antidotted samples 
the melting takes place in the zones with weaker pinning 
when compared to the case of not perforated samples in 
which the measured vortex phase transition takes place 
in the zones of intrinsic pinning. As a consequence of this 
also the hysteresis is much easier detectable in antidot- 
ted samples, in regions of the I — V characteristics not 
to close to the i?c2(T) curve. 

The influence of the regular array of antidots on the 
vortex properties is also confirmed by the behavior in 
our samples of the vortex correlation length in the liquid 
phase defined as |Q 

where 6 is a constant of the order of the unity, z/=0. 36963 
and f+o, being the smallest characteristic length scale in 
the liquid, is of the order of ao, the vortex lattice param- 
eter. According to the melting theory the shear viscosity 
ri{T) of the vortex liquid starts to grow up approacching 
the liquid-solid line from above, when r/ ~ @- 
The curves of ^+ versus the temperature for different 
applied magnetic fields are reported in figure || for the 
sample NCMH. All the curves start to diverge when the 
value of ^+ becomes comparable to the average distance 
among the antidots, i.e. when ^+ « d =1 /xm. Similar 
results have been obtained for all the samples investi- 
gated. This is exactly what we expect, if we look at the 
investigated system as a vortex ensemble constrained in 
the narrow channels among the lines of the antidots. In 
this case, infact, the melting transition at H = Hm is 
going to be observed when the correlation length of the 
liquid ^+ reaches the value of the width of the channels. 
The results shown in figure ^ clearly indicate the influ- 
ence of the antidots lattice on the vortex dynamics in our 
samples. 
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FIG. 8. Dependence of the correlation length ^+ versus 
the temperature at different magnetic fields for the sample 
NCMH. 

In conclusion, we have studied transport properties of 
superconducting (Nb)-spin glass (CuMn) multilayer with 
a regular array of antidots by measuring I — V curves in 
perpendicular magnetic fields. The measurements have 
been performed far above the matching conditions. The 
dynamic phase diagram has been drawn out from the 
analysis of these measurements. Two regions correspond- 
ing to plastic flux flow motion and to the motion of the 
vortex liquid have been distinguished. Melting occurs 
mostly due to quantum fluctuations and the presence of 
the antidots renders more easier to detect the melting 
due to the weaker pinning in the zones with higher local 
current density. 
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